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We discuss the “constant speed of sound” (CSS) parametrization of the equation of state of high-density mat¬ 
ter and its application to the field correlator method (FCM) model of quark matter. We show how observational 
constraints on the maximum mass and typical radius of neutron stars are expressed as constraints on the CSS 
parameters. We find that the observation of a 2M0 star already severely constrains the CSS parameters, and is 
particularly difficult to accommodate if the squared speed of sound in the high-density phase is assumed to be 
around 1/3 or less. 

We show that the FCM equation of state can be accurately represented by the CSS parametrization, which 
assumes a sharp transition to a high-density phase with density-independent speed of sound. We display the 
mapping between the FCM and CSS parameters, and see that FCM only allows equations of state in a restricted 
subspace of the CSS parameters. 

PACS numbers: 25.75.Nq, 26.60.-c, 97.60.Jd 


I. INTRODUCTION 

There are many models of matter at density significantly 
above nuclear saturation density, each with their own param¬ 
eters. In studying the equation of state (EoS) of matter in this 
regime it is therefore useful to have a general parametrization 
of the EoS which can be used as a generic language for relat¬ 
ing different models to each other and for expressing experi¬ 
mental constraints in model-independent terms. In this work 
we use the previously proposed “constant speed of sound” 
(CSS) parametrization KMl (for applications, see, e.g., El). 
We show how mass and radius observations can be expressed 
as constraints on the CSS parameters. Here we analyze a spe¬ 
cific example, where the high-density matter is quark matter 
described by a model based on the field correlator method 
(Sec. 0 , showing how its parameters can be mapped on to 
the CSS parameter space, and how it is constrained by cur¬ 
rently available observations of neutron stars. 

The CSS parametrization is applicable to high-density 
equations of state for which (a) there is a sharp interface be¬ 
tween nuclear matter and a high-density phase which we will 
call quark matter, even when (as in Sec. 0 we do not make 
any assumptions about its physical nature; and (b) the speed 
of sound in the high-density matter is pressure-independent 
for pressures ranging from the first-order transition pressure 
up to the maximum central pressure of neutron stars. One can 
then write the high-density EoS in terms of three parameters: 
the pressure ptrans of the transition, the discontinuity in energy 
density Ae at the transition, and the speed of sound Cqj^ in the 
high-density phase. Eor a given nuclear matter EoS eNM(p), 
the full CSS EoS is then 

£nm(p) P ^ Ptrans 

( ^NM (Ptrans ) T CQjyj {p Ptrans ) F ^ Ftrans 

( 1 ) 

The CSS form can be viewed as the lowest-order terms of 
a Taylor expansion of the high-density EoS about the tran¬ 
sition pressure. Eollowing Ref. fl], we express the three 


parameters in dimensionless form, as ptran.s/etrans, Ae/Ctrans 
(equal to A — 1 in the notation of Ref. ||5]) and Cgjyj, where 

£trans = £nM (Ptrans ) ■ 

The assumption of a sharp interface will be valid if, for ex¬ 
ample, there is a first-order phase transition between nuclear 
and quark matter, and the surface tension of the interface is 
high enough to ensure that the transition occurs at a sharp in¬ 
terface (Maxwell construction) not via a mixed phase (Gibbs 
construction). Given the uncertainties in the value of the sur¬ 
face tension fbHHl, this is a possible scenario. One can also 
formulate generic equations of state that model interfaces that 
are smeared out by mixing or percolation smi. 

The assumption of a density-independent speed of sound is 
valid for a large class of models of quark matter. The CSS 
parametrization is an almost exact fit to some Nambu-Jona- 
Lasinio models EIMII- The perturbative quark matter 
EoS lITSl also has roughly density-independent Cgj^, with a 
value around 0.2 to 0.3 (we use units where h = c = 1), above 
the transition from nuclear matter (see Eig. 9 of Ref. ifT^ ). 
In the quartic polynomial parametrization ini, varying the 
coefficient a 2 between ±(150MeV)^, and the coefficient aq 
between 0.6 and 1, and keeping ntians/«o above 1.5 (mq = 
0.16 fm^^ is the nuclear saturation density), one finds that Cg|^ 
is always between 0.3 and 0.36. It is noticeable that mod¬ 
els based on relativistic quarks tend to have Cgjyj « close to 
1/3, which is the value for systems with conformal symmetry, 
and it has been conjectured that there is a fundamental bound 
CgM <1/3 lITSl . although some models violate that bound, 
e-g- IISlISOl or l[T4l (parametrized in || 2 l). 

In Sec. [I^ we show how the CSS parametrization is con¬ 
strained by observables such as the maximum mass the 

radius of a maximum-mass star, and the radius Ri 4 of a star 
of mass 1 . 4 M 0 . In Secs. [In] - [IV| we describe a specific model, 
based on a Brueckner-Hartree-Eock (BHE) calculation of the 
nuclear matter EoS and the field correlator method (ECM) for 
the quark matter EoS. We show how the parameters of this 
model map on to part of the CSS parameter space, and how 
the observational constraints apply to the ECM model param- 
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eters. Section[V]gives our conclusions. 


II. CONSTRAINING THE CSS PARAMETERS 
A. Topology of the mass-radius relation 



FIG. 1: (Color online). Schematic phase diagram (from UJ) for hy¬ 
brid star branches in the mass-radius relation of compact stars. We 
fix and vary ptrans/Ctrans and Ae/ enans- The four regions are (A) 

no hybrid branch (“absent”); (B) both connected and disconnected 
hybrid branches; (C) connected hybrid branch only; and (D) discon¬ 
nected hybrid branch only. 

We use the term “hybrid star” to refer to stars whose cen¬ 
tral pressure is above ptrans. and so they contain a core of 
the high-density phase. The part of the mass-radius relation 
that arises from such stars is the hybrid branch. In all mod¬ 
els of nuclear/quark matter we find the same four topologies 
of the mass-radius curve for compact stars; the hybrid branch 
may be connected to the nuclear branch (C), or disconnected 
(D), or both may be present (B) or neither (A). The occur¬ 
rence of these as a function of the CSS parameters ptrans/etrans 
and Ae/Ctrans at fixed Cgjyj is shown schematically in Fig. 
(taken from Ref. m). The mass-radius curve in each region 
is depicted in inset plots, in which the thick green line is the 
hadronic branch, the thin solid red lines are stable hybrid stars, 
and the thin dashed red lines are unstable hybrid stars. 

In the phase diagram the solid red line shows the thresh¬ 
old value ACcrit below which there is always a stable hybrid 
star branch connected to the neutron star branch. This critical 
value is given by ||5] [21] l22l 

■^£crit _ 1 _j_ 3 Ptrans ^ 2 ^ 

Grans ^ ^ Grans 

and was obtained by performing an expansion in powers of 
the size of the core of high-density phase. Equation (|^ is 
an analytic result, independent of Cqj^ and the nuclear matter 
EoS. The dashed and dot-dashed black lines mark the border 
of regions where the disconnected hybrid star branch exists. 


The position of these lines depends on the value of Cqj^ and 
(weakly) on the accompanying nuclear matter EoS H]. 

Once a nuclear matter EoS has been chosen, any 
high-density EoS that is well approximated by the CSS 
parametrization can be summarized by giving the values of the 
three CSS parameters, corresponding to a point in the phase 
diagram. We then know what sort of hybrid branches will be 
present. 


B. Maximum mass of hybrid stars 


Property 

BHF, Avig 
- 1 - UVIX TBF 

DBHF, 
Bonn A 

Saturation baryon density no(fm^’^) 

0.16 

0.18 

Binding energy/baryon E/A (MeV) 

-15.98 

-16.15 

Compressibility Kq (MeV) 

212.4 

230 

Symmetry energy Sq (MeV) 

31.9 

34.4 

L = 3«o [dSa/dnU (MeV) 

52.9 

69.4 

Maximum mass of star (M©) 

2.03 

2.31 

Radius of the heaviest star (km) 

9.92 

11.26 

Radius of M = I. 4 M 0 star (km) 

11.77 

13.41 


TABLE I: (Color online). Calculated properties of symmetric nu¬ 
clear matter for the BHF and Dirac-Brueckner-Hatree-Fock (DBHF) 
nuclear equations of state used here. BHF is softer, and DBHF is 
stiffer (see Sec. |ng 

In Eig. 1^ we show how mass measurements of neutron 
stars can be expressed as constraints on the CSS parameters. 
Each panel shows dependence on ptrans/Grans and Ae/eo-ans 
for fixed Cqj^j, as in Eig.j^ The region in which the transition 
to quark matter would occur below nuclear saturation density 
(«tians < no) is excluded (hatched band at the left end) because 
in that region bulk nuclear matter would be metastable. There 
is also an upper limit on the transition pressure, which is the 
central pressure of the heaviest stable nuclear matter star. This 
depends on the hadronic EoS that had been assumed. 

The contours show the maximum mass of a hybrid star 
as a function of the EoS parameters. The region inside the 
M = 2Mq contour corresponds to EoSes for which the max¬ 
imum mass is less than 2 M 0 so it is shaded to signify that 
this region of parameter space for the high-density EoS is ex¬ 
cluded by the observation of a star with mass 2M© ll23l . Eor 
high-density EoSs with Cgj^ = 1 (right-hand plots), this re¬ 
gion is not too large, and leaves a good range of transition 
pressures and energy density discontinuities that are compat¬ 
ible with the observation. However, for high-density matter 
with Cgj^ = 1/3 (left-hand plots), which is the typical value in 
many models (See Sec.|^, the M^ax > 2 M 0 constraint elim¬ 
inates a large region of the CSS parameter space lEIIl. We 
discuss this in more detail below. 

The upper plots in Eig.j^are for a stiffer nuclear matter EoS, 
dbheEtI. and the lower plots are for a softer nuclear matter 
EoS, BHE ||25]| (see Sec. [In|). Properties of these nuclear mat¬ 
ter EoSs are given in Table Ifi As one would expect, the stiffer 
EoS gives rise to heavier (and larger) stars, and therefore al- 
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FIG. 2: (Color online). Contour plots showing the maximum hybrid star mass as a function of the CSS parameters of the high-density EoS. 
Each panel shows the dependence on the CSS parameters ptrans/Cti-ans and Ae/eti-ans- The left plots are for Cgiyj = 1/3, and the right plots are 
for Cgjyj = 1. The top row is for a DHBE (stiff) nuclear matter EoS, and the bottom row is for a BHE (soft) nuclear matter EoS. The grey 
shaded region is excluded by the measurement of a 2M0 star. The hatched band at low density (where utrans < no) is excluded because bulk 
nuclear matter would be metastable. The hatched band at high density is excluded because the transition pressure is above the central pressure 
of the heaviest stable hadronic star. 


lows a wider range of CSS parameters to be compatible with 
the 2 M 0 measurement. 

In Fig. 1^ the dot-dashed (red) contours are for hybrid stars 
on a connected branch, while the dashed (blue) contours are 
for disconnected branches. As discussed in Ref. m, when 
crossing the near-horizontal boundary from region C to B 
the connected hybrid branch splits into a smaller connected 
branch and a disconnected branch, so the maximum mass of 
the connected branch smoothly becomes the maximum mass 
of the disconnected branch. Therefore the red contour in the 
C region smoothly becomes a blue contour in the B and D 
regions. When crossing the near-vertical boundary from re¬ 
gion C to B a new disconnected branch forms, so the con¬ 
nected branch (red dot-dashed) contour crosses this boundary 


smoothly. 

In each panel of Fig. the physically relevant allowed re¬ 
gion is the white unshaded region. The grey shaded region is 
excluded by the existence of a 2 M 0 star. We see that increas¬ 
ing the stiffness of the hadronic EoS or of the quark matter 
EoS (by increasing Cqj^) shrinks the excluded region. 

Eor both the hadronic EoSs that we study, the CSS param¬ 
eters are signihcantly constrained. Erom the two left panels 
of Eig. 1^ one can see that if, as predicted by many models, 
Cqm ^ 1/3, then we are limited to two regions of parame¬ 
ter space, corresponding to a lowpressure transition or a high 
pressure transition. In the low-transition-pressure region the 
transition occurs at a fairly low density ntrans ^ 2 no, and a 
connected hybrid branch is possible. In the high-transition- 
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pressure region the connected branch (red dot-dashed) con¬ 
tours are, except at very low Ae, almost vertical, correspond¬ 
ing to EoSs that give rise to a very small connected hybrid 
branch which exists in a very small range of central pressures 
Pcent just above ptians- The maximum mass on this branch is 
therefore very close to the mass of the purely hadronic matter 
star with pcent = ftians- The mass of such a purely hadronic 
star is naturally independent of parameters that only affect the 
quark matter EoS, such as Ae and CQ]y[, so the contour is ver¬ 
tical. These hybrid stars have a tiny core of the high-density 
phase and cover a tiny range of masses, of order 10 ^^ M© or 
less, and so would be very rare. 

Disconnected hybrid branches are of special interest, be¬ 
cause they give a characteristic signature in mass-radius mea¬ 
surements. Eor both the hadronic EoSs that we study, the re¬ 
gion B and D, where disconnected hybrid star branches can 
occur, are excluded for Cqj^ ^1/3. Even for larger Cgiyj dis¬ 
connected branches only arise if the nuclear matter EoS is suf¬ 
ficiently stiff. It is interesting to note that using an extremely 
stiff hadronic matter EoS such as DD2-EV ll26l can further 
shrink the region that is excluded by the M^ax > 2 Mq con¬ 
straint, allowing disconnected branches of hybrid stars to oc¬ 
cur. 


C. Minimum radius of hybrid stars 

In Eig. we show contour plots of the radius of the 
maximum-mass star (on either a connected or disconnected 
hybrid branch) as a function of the CSS quark matter EoS pa¬ 
rameters. Since the smallest hybrid star is typically the heav¬ 
iest one, this allows us to infer the smallest radius that arises 
from a given EoS. 

The layout is as in Eig. each panel shows dependence on 
Ptrans/etrans and Ag/gtrans for fixed CQ]y.j; the plots on the left 
are for Cgjyj = 1/3 and the plots on the right are for Cgp^ = 1; 
the plots on the top are for the stiffer DBHE nuclear matter 
EoS, while the lower plots are for the softer BHE nuclear mat¬ 
ter EoS. As in Eig.|^ the region that is eliminated by the ob¬ 
servation of a 2 M 0 star is shaded in grey. 

The smallest stars, with radii as small as 9 km, occur when 
the high-density phase has the largest possible speed of sound 
Cqm = 1. They are disconnected branch stars arising from 
EoSs having a low transition pressure (utrans ^ 2no) with a 
fairly large energy density discontinuity (Ag/go-ans ^ !)■ 

As in Eig. the contours in the high-transition-pressure 
region are almost vertical because the hybrid branch is then a 
very short extension to the nuclear mass-radius relation, and 
its radius is close to that of the heaviest purely hadronic star, 
which is independent of Ag/gtrans Cgj^. The radius of the 
hybrid stars decreases with ptrans in this region, because the 
radius of hadronic stars decreases with central pressure. 

Eor Cgjy^j = 1 /3, the allowed low-transition-pressure region 
is disconnected from the high-transition-pressure region and 
is so small that it is hard to see on this plot. By magnifying 
it (left-hand plots of Fig. we see that in this region the ra¬ 
dius contours closely track the border of the allowed region 


(the Mmax = 2 M 0 line) so we can say that the radius must 
be greater than 11.5 km almost independent of the transition 
pressure and hadronic EoS. For a stiff hadronic EoS this mini¬ 
mum is raised to 11.7 km. These values are comparable to the 
minimum radius of about 11.8 km found in Ref. ifTSl . which 
explored a larger set of hadronic EoSs but did not explore the 
full CSS parameter space for the high-density EoS. If a star 
with radius smaller than this minimum value were to be ob¬ 
served, we would have to conclude that either the transition 
occurs outside the low-density region or that Cgjyj is greater 
than 1/3. In the magnified figure we also show how the ex¬ 
cluded region would grow if a 2 . 1 M 0 star were to be observed 
(long-dashed line for connected branch stars and short-dashed 
line for disconnected branch stars). This would increase the 
minimum radius to about 12.1 km for the soft hadronic EoS 
and 12.2km for the stiff hadronic EoS. 


D. Typical radius of hybrid stars 

In Fig. 1^ we show contours (the U-shaped lines) of typical 
radius of a hybrid star, defined as 4 , the radius of a star of 
mass I. 4 M 0 , as a function of the CSS quark matter EoS pa¬ 
rameters. The contours only fill the part of the CSS parameter 
space where there are hybrid stars with that mass. The dashed 
(magenta) lines delimit that region which extends only up to 
moderate transition pressure. 

The overall behavior is that, at fixed Ag/gtrans. the typi¬ 
cal radius is large when the transition density is at its low¬ 
est. As the transition density rises the radius of a I. 4 M 0 
star decreases at first, but then increases again. This is re¬ 
lated to the previously noted fact that when one fixes the 
speed of sound of quark matter and increases the bag con¬ 
stant (which increases ptrans/Ctrans and also varies Ag/gnans in 
a correlated way) the resultant family of mass-radius curves 
all pass through the same small region in the M-R plane: the 
M{R) curves “rotate” counterclockwise around this hub (see 
Fig. 2 of Ref. 1221). In our case we are varying ptrans/fitrans at 
fixed Ag/gtrans, SO the hub itself also moves. At low transi¬ 
tion density the hub is below 1.4M0, so R\ a, decreases with 
Ptrans/etrans- At high transition density the hub is at a mass 
above I. 4 M 0 so Ri ,4 will increase with ptrans/Ctrans- 

The smallest stars occur for Cgjyj = 1 (right-hand plots), 
where 4 > 9.5km at large values of the energy density 
discontinuity, and the radius rises as the discontinuity is de¬ 
creased. This is consistent with the absolute lower bound of 
about 8.5 km ESl for the maximally compact Cgj^j = 1 star 
obeying Mmax > 2 M 0 . 

For Cgp^ = 1/3 the allowed region at low transition pressure 
is small, so in the right panels of Fig. we show a magnifica¬ 
tion of this region. We see that in the allowed (Mmax > 2 M 0 
and Htrans > Ho) region there is a minimum radius 12 . 2 km 
for the BHE (soft) hadronic EoS, and about 12.5km for the 
DBHE (stiff) hadronic EoS. This minimum is attained at the 
lowest possible transition density, Wtrans ~ no- As the transi¬ 
tion density rises to values around 2 no, the minimum radius 
rises to 12.5 (BHE) or 13.3 km (DBHE). This is comparable 
to the minimum radius of about 13 km found in Ref. m. 
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FIG. 3: (Color online). Contour plots showing the radius of the maximum-mass star as a function of the CSS parameters. Dashed lines are for 
the case where this star is on the disconnected branch; for dot-dashed lines it is on the connected branch. The grey shaded region is excluded 
by the measurement of a 2Mq star. The hatched band at low density (where «trans < no) is excluded because bulk nuclear matter would be 
metastable. The hatched band at high density is excluded because the transition pressure is above the central pressure of the heaviest stable 
hadronic star. For a magnified version of the low-transition-pressure region for Cq^ = 1/3, see Fig.l^ 


which explored a wider range of hadronic EoSs but assumed 
Htrans = 2no. These results are consistent with the lower bound 
on Ri 4 for Cqj^ = 1/3 of about 11 km established in Ref. Il28l 
(Fig. 5) using the EoS that yields maximally compact stars 
(corresponding to CSS with ptrans = 0 and Cgjy^j = 1 /3) obey¬ 
ing Mmax > 2 M 0 . If a 1 . 4 M 0 star were observed to have 
radius below the minimum value, one would have to conclude 
that either it is not a hybrid star or that Cgjyj >1/3. 


III. THE BHF AND DBHF EOS OF NUCLEAR MATTER 


We now discuss in more detail the nuclear matter equa¬ 
tions of state that we use in this work. We adopt the BHF 
scheme, in which the only input needed is the realistic free 
nucleon-nucleon (NN) interaction V in the Brueckner-Bethe- 
Goldstone (BBG) equation for the reaction matrix G, 


G[p;td] =y-F Y, ^ 


(O - e{ka) - e{kb) 


G[p-,(o], 


(3) 


The dashed line shows how the excluded region would grow 
if a star of mass 2.1 M 0 were to be observed. This would 
increase the minimum radius to about 12.7 (BHF) or 13 km 
(DBHF). 


where p is the nucleon number density, and co the start¬ 
ing energy. The propagation of intermediate baryon pairs 
is determined by the single-particle energy e(k\p) = ^ + 
U{k-,p), and the Pauli operator Q. Because of the oc- 
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FIG. 4: (Color online). Contour plots similar to Fig. showing the radius of a hybrid star of mass M = 1.4M0 as a function of the CSS 
parameters. Such stars only exist in a limited region of the space of EoSs [delimited by dashed (magenta) lines]. The grey shaded region 
is excluded by the observational constraint Mmax > 2Mq. For a magnified version of the low-transition-pressure region for Cgjyj = 1/3, see 
Fig.[5l 


currence of G in the single-particle potential U{k',p) = 
Re^j.'</if {kk'\G[p',e{k)+e{k!)]\kk!)a, where the subscript “a” 
indicates antisymmetrization of the matrix element, the BBG 
equation (Eq. Q) has to be solved in a self-consistent manner 
for several momenta of the particles involved, at the consid¬ 
ered densities. 

In the nonrelativistic BHF approximation the energy per nu¬ 
cleon is given by ll29l 

The nuclear EoS can be calculated with good accuracy in the 
Brueckner two hole-line approximation with the continuous 
choice for the single-particle potential, and the results in this 
scheme are quite close to the calculations which include also 
the three hole-line contribution Il30l432l . The dependence on 


the NN interaction, also within other many-body approaches, 
has been systematically investigated in Ref. Il3?l . 

It is well known that, in order to reproduce the correct sat¬ 
uration point of symmetric nuclear matter, we must introduce 
nuclear three-body forces (TBFs). In our approach the TBF is 
reduced to a density-dependent two-body force by averaging 
over the position of the third particle, assuming that the prob¬ 
ability of having two particles at a given distance is reduced 
according to the two-body correlation function 011351. 

In this work we use the Argonne vig NN potential 1^ . and 
the so-called Urbana model for TBFs, which consists of an at¬ 
tractive term due to two-pion exchange with excitation of an 
intermediate A resonance, and a repulsive phenomenological 
central term OTl OBJ. Those TBFs produce a shift of about 
+ 1 MeV in energy and —0.01 fm in density. This adjust¬ 
ment is obtained by tuning the two parameters contained in the 
TBFs, and was performed to get an optimal saturation point 
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1.4M0 star. The region under and to the left of the hatched bar is probably unphysical because ntrans < and it was excluded (hatched band) 
in earlier figures. The grey shaded region is excluded by the observational constraint A/max > TMq. The dashed line shows how that region 
would grow if a 2.1 Mq star were observed. 


Il3l|35l. At present the theoretical status of microscopically 
derived TBFs is still quite rudimentary; however a tentative 
approach has been proposed using the same meson-exchange 
parameters as the underlying NN potential Il^l40l . 

Along with the nonrelativistic BHF EoS we consider its rel¬ 
ativistic counterpart, the DBHF scheme ll24l where the Bonn- 
A potential is used for the nucleon-nucleon interaction. In 
the low density region (p < 0.3 fm^^), the BHF (including 
TBF) and DBHF equations of state are very similar, whereas 
at higher densities the DBHF is slightly stiffen The discrep¬ 
ancy between the nonrelativistic and relativistic calculation 
can be easily understood by recalling that the DBHF treatment 
is equivalent to introducing in the nonrelativistic BHF the TBF 
corresponding to the excitation of a nucleon-antinucleon pair, 
the so-called Z-diagram ED, which is repulsive at all densi¬ 
ties. In the BHF treatment with Urbana TBF, both attractive 


and repulsive TBF are introduced and therefore a softer EoS 
is expected. We report in Table |I] the main properties of both 
EoSs. 

In this work we perform all calculations for both the BHE 
and DBHE equations of state for hadronic matter. This pro¬ 
vides a reasonable range of possible hadronic EoSs, and al¬ 
lows us to gauge the sensitivity of our results to this source of 
uncertainty, although even DBHE is not as stiff as an ultrastiff 
hadronic EoS such as DD2-EV ll26ll (see Sec. |IIB| l. 

We do not include the effects of hyperons because these are 
unknown, and including them would not increase the physical 
accuracy of our results. Calculating a hyperonic EoS requires 
knowledge of hyperon interactions with other baryons, and 
there is little data on hyperon-nucleon interactions ll42l and 
none on hyperon-hyperon interactions or three-body interac¬ 
tions. There have been various conjectures about the hyperon 
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interaction ll4?H46l and how to include it in BHF ll44ll47ll48l baryon density ps ll55l . in nuclear matter. However, for sim- 
and DBHF ll49l but there is no consensus on the correct result. plicity, in the following we treat both Vi and G 2 as numerical 

parameters with no dependence on pg. 


IV. QUARK MATTER VIA THE EIELD CORRELATOR 
METHOD 


B. The ECM EoS and the CSS parametrization 


A. The ECM EoS 

The approach based on the FCM provides a natural treat¬ 
ment of the dynamics of confinement in terms of the color 
electric {D^ and Df) and color magnetic {D^ and D^) Gaus¬ 
sian correlators, the former being directly related to confine¬ 
ment, so that its vanishing above the critical temperature im¬ 
plies deconfinement Eq). The extension of the FCM to fi¬ 
nite temperature T at chemical potential Pq = 0 gives analyt¬ 
ical results in reasonable agreement with lattice data, giving 
us some confidence that it correctly describes the deconfine¬ 
ment phase transition ifSTl l5^ . In order to derive an EoS of 
the quark-gluon matter in the range of baryon density typical 
of the neutron star interiors, we have to extend the FCM to 
nonzero chemical potential iniEii. In this case, the quark 
pressure for a single flavor is simply given by 






(5) 


where 


(/)v(a) = J du [u ^/^exp 


\J iP- — a -fl^ 

( 6 ) 


with V = niq/T , and Vi is the large distance static qq poten¬ 
tial whose value at zero chemical potential and temperature is 
V\{T = Pb—0) = 0.8 to 0.9GeV |l53][54l. The gluon contri¬ 
bution to the pressure is 



and the total pressure is 


(7) 


E P^+Pg 

j=U,d,S 


{n-lNf) G2 
32 2 


( 8 ) 


where and Pg are given in Eqs. ^ and Q, and Nf 
is the number of flavors. The last term in Eq. corre¬ 
sponds to the difference of the vacuum energy density in the 
two phases, G 2 being the gluon condensate whose numeri¬ 
cal value, determined by the QCD sum rules at zero tempera¬ 
ture and chemical potential, is known with large uncertainty, 
G 2 = 0.012 ±0.006 GeV^. At finite temperature and vanish¬ 
ing baryon density, a comparison with the recent available lat¬ 
tice calculations provides clear indications about the specific 
values of these two parameters, and in particular their values 
at the critical temperature Tc- Some lattice simulations suggest 
no dependence of Vi on p^, at least for very small pg, while 
different analyses suggest a linear decreasing of G 2 with the 


The CSS parametrization will be applicable to the ECM 
EoS if the speed of sound in the ECM EoS depends only 
weakly on the density or pressure. In Eig. we show that 
this is indeed the case. The upper panel shows the speed of 
sound vs. pressure in the ECM quark matter EoS for different 
values of the ECM parameters, displayed in the lower panel. 
We see that the speed of sound varies by less than 5% over 
the considered range of pressures along each curve, and lies 
in the interval 0.28 < CQjyj <1/3. The value of Cgjyj shows a 
weak dependence on Vi and extremely weak dependence on 
G 2 , which appears as an additive constant in the quark mat¬ 
ter EoS according to Eq. (§. The transition pressure is more 
sensitive to the ECM parameters, increasing rapidly with Vi 
and with G 2 . The energy density at a given pressure increases 
slightly with an increase in Vi or G 2 . 



FIG. 6: (Color online). The squared speed of sound Cqj^ [panel (a)] 
is displayed vs quark matter pressure for several values of Vi (in 
MeV) and G 2 (in GeV^). In panel (b), the FCM energy density is 
represented by full symbols, whereas the full lines denote the CSS 
parametrization given by Eq. Q. 

To illustrate how well the CSS parametrization fits the FCM 
EoS, we show in the lower panel of Fig. that, for the same 
FCM parameter choices, we can always And suitable values of 
the CSS parameters which fit the FCM calculation extremely 
well. This means that there exists a well-defined mapping be¬ 
tween the FCM parameters (Vi, G 2 ) and the CSS parameters 
(Ptrans/etians, Ac^trans, ^ote that the mapping depends 

on the EoS of the hadronic matter. 

The mapping is displayed in Fig.|^ which shows the region 
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BHF (soft) NM, FCM 


DBHF (stiff) NM, FCM 



FIG. 7: (Color online). The mapping of the FCM quark matter model onto the CSS parametrization. Results are obtained using the BHF 
(left panel) and DBHF (right panel) nuclear matter EoS. The undecorated curves are the phase boundaries for the occurrence of connected and 
disconnected hybrid branches (compare Figs.[T]an d[g. The thin dashed (black) line and the solid (black) line studded with circles delimit the 
region yielded by the FCM model. Within that region, lines decorated with symbols give CSS parameter values for FCM quark matter as G 2 
is varied at constant Vi (given in MeV). The (orange) cross denotes the EoS with the highest ptrans, which gives the heaviest FCM hybrid star. 
See the text for details. 


of the CSS parameter space where FCM equations of state 
are found. As in the phase diagrams in Sec. we show the 
plane whose coordinates are the CSS parameters Ae/Ctrans and 
Ptrans/Ctrans- For the hadronic EoS we use BHF (left panel) and 
DBHF (right panel). The lines without points represent the 
phase boundaries, as for the figures in Sec.|^ for connected 
and disconnected branches. Whether a given FCM EoS yields 
stable hybrid stars depends on which of those phase regions 
(see Fig. it is in. The solid (green) phase boundary with a 
cusp at ptran.s/etrans ~ 0.17 delimits the region with a discon¬ 
nected branch for Cgjyj = 1 /3, while the nearby dashed (green) 
line is for CQjyj = 0.28, so these span the range of Cgjyj relevant 
for the FCM, as discussed in Fig. It is evident that the de¬ 
pendence on Cgjyi is tiny and negligible for practical purposes. 

The thin dashed (black) line and the solid (black) line stud¬ 


ded with circles delimit the equations of state yielded by the 
FCM calculation. Within that region, the lines studded with 
points show the CSS parametrization of the FCM quark mat¬ 
ter EoS, where along each line we keep Vi constant and vary 
G 2 . Above that region, which corresponds to negative values 
of Vi, the EoS cannot reproduce the 2 M 0 limit, and in this 
sense is unphysical (Fig. [^. Below that region, there would 
be no transition from hadronic to quark matter, as explained 
below. 

In Fig. [7] Vi varies from 0 up to the maximum value at 
which hybrid star configurations occur, which is indicated 
by an (orange) cross. For the BHF case that value is Vi = 
240 MeV, G 2 = 0.0024 GeV'^ and for the DBHF case it is 
Vi = 255 MeV, G 2 = 0.0019 GeV'^. Along each FCM curve in 
Fig.|^the parameter G 2 starts at the minimum value at which 
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BHF(soft) NM, FCM 


1.0 2.0 3.0 4.0 ° 6.0 7.0 



DBHF (stiff) NM, FCM 


n, /n„ 

trans 



FIG. 8: (Color online). Contour plots, analogous to Fig.[^ showing the maximum mass of hybrid stars with FCM quark matter cores, given in 
terms of the corresponding CSS parameter values rather than the original FCM parameter values. As in Fig.|^ solid lines are phase boundaries 
(compare Figs. [T]and|^. The shaded sectors indicate the parameter regions accessible by the FCM and with Mmax > TMq. In each panel the 
lower border of the shaded region meets the phase boundary (red line) at the point with highest value of Vi reported in Fig. as an orange 
cross. Note the different scales on the x axis for the two panels. 


there is a phase transition from hadronic to FCM quark mat¬ 
ter; at lower G 2 the quark and the hadronic pressures p(/r) do 
not cross at any /r. On each curve one point is labeled with its 
value of 02 /( 10 ^^ GeV^), and subsequent points are at inter¬ 
vals where G 2 increases in increments of 1 in the same units. 

We observe that along each line of constant Vi, ptrans/fitrans 
grows with G 2 . This can be explained by recalling the linear 
dependence of the quark pressure on G 2 in Eq. (|^, so that, at 
fixed chemical potential, an increase of G 2 lowers the quark 
pressure, making quark matter less favorable, and shifting the 
transition point to higher chemical potential or pressure. This 
was already discussed in Ref. Il56l for BHF nuclear matter, 
and is equally applicable to DBHF nuclear matter. Obviously 
if G 2 becomes too small, the phase transition takes place in a 
region of low densities where finite nuclei are present, and the 
homogeneous nuclear matter approach becomes invalid. 


The qualitative behavior of the curves of constant Vi can 
be understood in terms of the Maxwell construction between 
the purely hadronic phase and the quark phase. The fact 
that Ae/Ctrans goes through a minimum (which is always at 
Ptrans/Ctrans « 0.1) as G 2 is increased at constant Vi can be un¬ 
derstood from Fig. 2 of Ref. ESI, which shows pressure p as 
a function of baryon density n and the location of the hadron 
(BHF EoS) to quark (FCM EoS) transition when G 2 is varied. 
The hadronic EoS is strongly curved, especially at low pres¬ 
sure, while the FCM EoS is closer to a straight line. Conse¬ 
quently, the baryon density difference between the two phases 
at a given pressure has a minimum at densities around 2 no, 
which corresponds to ptrans/Ctrans ~ 0.1. As G 2 increases, the 
transition pressure rises, scanning through this minimum. It 
follows that the energy density difference also goes through a 
minimum, because e = fin — p, and p and fi are continuous at 
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FIG. 9: (Color online). Shaded areas show range of radii of stars 
with a given maximum mass when varying FCM parameters. The 
thick dashed lines indicate the purely hadronic mass-radius configu¬ 
rations. The crosses at the top of the shaded regions correspond to 
the maximal configuration indicated by the same symbol in Fig. 
The observational constraint da on the star mass is indicated by a 
horizontal line. 


the transition, so Ae = /r An. The DBHF hadronic EoS is very 
similar to BHF at low pressure, so the curves have their min¬ 
ima at the same value of ptrans/Ctrans in both panels of Fig.|7] 
We also see in Fig.|7]that an increase of Vi moves the curves 
slightly downward and to the right. This is expected since Vi 
is a measure of the interparticle strength, and is inversely pro¬ 
portional to the pressure of the system, so that the pressure 
decreases as Vi is increased at fixed /r, and, as already dis¬ 
cussed for the parameter G 2 , a decrease of the quark pressure 
raises ptrans- The role of Vi and Gi in the quark EoS discussed 
so far, provides in the same way a qualitative understanding 
of CQ]y[ in panel (a) in Fig. although, as already noticed, the 
effect in Fig.j^of the change in Cqj^j is negligible. 


C. Expected properties of mass-radius curves 

By comparing Fig. |7] with Fig. [T] we can see that when 
combining FCM quark matter with BHF (soft) nuclear mat¬ 
ter, the physically allowed range of FCM parameter values 
yields EoSs that are mostly in regions C and A, where there is 
no disconnected hybrid branch. At the lowest transition den¬ 
sities the FCM EoS can achieve a large enough energy density 
discontinuity to yield a disconnected branch (region D). 

For the DBHF (stiff) nuclear EoS there is a wider range 
of values of Vi and Gj that give disconnected branches, and 


some of them give simultaneous connected and disconnected 
branches. This difference can be understood in terms of the 
stiffness of the EoSs. A change from a soft hadronic EoS 
(BHF) to a stiff one (DBHF) produces a steeper growth of the 
hadronic pressure as a function of the baryon density. Re¬ 
ferring again to Fig. 2 of Ref fSb), this pulls the DBHF p{n) 
curve further away from the FCM curve, giving a larger differ¬ 
ence in baryon density at a given pressure, and hence, as noted 
above, a larger Ae. This is why the curves for DBHFh-FCM 
(right panel of Fig. are shifted upwards along the Ae/euans 
axis compared to the BHFh-FCM curves (left panel of Fig. 0. 

We can calculate the maximum mass of a hybrid star con¬ 
taining a FCM core as a function of the FCM parameters, and 
then use the mapping described above to obtain the CSS pa¬ 
rameter values for each FCM EoS, producing a contour plot of 
maximum mass (Fig. [8|l for BHF (left panel) and DBHF (right 
panel) hadronic EoS. Given that the CSS parametrization is a 
fairly accurate representation of the FCM EoS, one would ex¬ 
pect this to be very similar to the corresponding plot for CSS 
itself with Cgj^ = 1/3 (Fig |^i and this is indeed the case. The 
contours in Fig.|^are restricted to the region corresponding to 
physically allowed FCM parameter values, so they end at the 
edges of that region. 

The triangular shaded area at the edge of each panel shows 
the region of the parameter space that is accessible by the 
FCM and is consistent with the measurement of a 2 M 0 , by 
having hybrid stars of maximum mass greater than 2M©. 
The (orange) cross in each panel of Fig. is at the high- 
transition-pressure corner of that triangular area. The heav¬ 
iest BHFh-FCM hybrid star has a mass of 2.03M©, and the 
heaviest DBHFh-FCM hybrid star has a mass of 2.3IM©. 

As noted in Sec. Inl the hybrid stars in this physically al¬ 
lowed and FCM-compatible region of the phase diagram lie 
on a very tiny connected branch, covering a very small range 
of central pressures and masses and radii, and would therefore 
occur only rarely in nature. These stars have very small quark 
matter cores (see Ref. 11, Figs. 5-6), and their mass and ra¬ 
dius are very similar to those of the heaviest purely hadronic 
star, but there could be other clear signatures of the presence 
of the quark matter core, such as different cooling behavior. 

The CSS parametrization has another region where heavy 
hybrid stars occur, at low transition pressure (see Fig. 1^, but 
the FCM does not predict that the quark matter EoS could be 
in that region. 

To characterize the radius of FCM hybrid stars we can¬ 
not construct contour plots like Fig. because, as we have 
just seen, the FCM predicts that only hybrid stars with mass 
very close to the maximum mass are allowed. There are no 
FCM hybrid stars with mass around 1.4M©. Instead, in Fig.[^ 
we show the range of radii of stars with a given maximum 
mass when varying FCM parameters, for our two different 
hadronic EoSs. The right-hand edge of each shaded region 
traces out the mass-radius relation for hadronic stars with the 
corresponding hadronic EoS. The FCM hybrid stars form very 
small connected branches which connect to the nuclear mat¬ 
ter where the central pressure reaches the transition pressure 
(see Sec. IIB 1 , so the hybrid stars do not deviate very far from 
the hadronic mass-radius curve. Hence the shaded regions in 
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Fig. 1^ are narrow, especially in the observationally allowed 
(iWmax > 2 M 0 ) region, which perfectly matches the predic¬ 
tion of CSS parametrization on the maximum-mass star ra¬ 
dius in the high-transition-pressure region (see the left panels 
of Fig.[^. For BHF (soft) nuclear matter, the hadronic stars, 
and hence the hybrid stars, are smaller because the nuclear 
mantle is more compressed by the self gravitation of the star. 

V. CONCLUSIONS 

We have shown how observational constraints on the mass 
and radius of hybrid stars can be expressed as constraints on 
the parameters of the CSS parametrization of the high-density 
EoS, which, in the space of possible models of quark matter, is 
a reasonably model-independent parametrization. Of course, 
physical predictions from CSS depend on the hadronic EoS 
with which it is combined. The CSS parametrization assumes 
a sharp transition from nuclear matter to a high-density phase 
such as quark matter, and that the speed of sound in that phase 
is independent of the pressure. We found that the observation 
of a 2Mq star constrains the CSS parameters significantly IT] 

M- 

If, as predicted by many physical models of quark matter, 
Cqm ^1/3, then for typical models of hadronic matter such 
as BHF or DBHF there are two possible scenarios (we discuss 
ultrastiff hadronic EoSs below). 

First, there is a low-transition-pressure scenario, where the 
transition to the high-density phase occurs at Utrans ^ 2no (the 
unshaded region on the left side of the two left panels of 
Figs. |2jQ. In this scenario, the hybrid branch of the mass- 
radius relation will be connected to the nuclear branch. In 
the Cgjyj ^1/3 and low-transition-pressure scenario there are 
strong constraints on the radius of the star, as shown in Fig. 
The radius of the maximum-mass star (which is typically the 
smallest possible star) must be greater than about 11.5 km, and 
the radius of a 1. 4 M 0 star must be greater than about 12.2 km 
1281 . For a stiffer hadronic EoS, these minima are raised by 
about 0.15 to 0.3 km. If a neutron star of mass 2 .IM 0 were 
observed then this constraint would tighten, increasing the 
minimum radius to about 12.1 km. If a star smaller than the 
minimum radius were observed, we would have to conclude 
that either the transition is outside the low-density regime 
or CQ]y[ > 1/3. Conversely, if theoretical considerations es¬ 
tablished that Cqj^ is smaller than 1/3, the minimum radius 
would become larger im. 

Secondly, there is a high-transition-pressure scenario (the 
white region on the right side of the left panels of Figs. m- 
This tends to give a very small branch of hybrid stars with 
tiny quark matter cores, occurring in a narrow range of cen¬ 
tral pressures just above the transition pressure. This is why 
the mass and radius contours become almost vertical in this 
region: the hybrid star has almost the same mass and radius 
as the heaviest purely hadronic star (the one where the central 
pressure is ptrans), and so the properties of these hybrid stars 
depend on the hadronic EoS (see Eig. 1^ via Ptrans/^^trans btlt 
not on quark matter properties such as Ae or Cgjyj. 

If the hadronic matter is extremely stiff (e.g. DD2-EV ll26l l 


or the quark matter has Cgjyj larger than 1/3 then a larger re¬ 
gion of the CSS parameter space becomes allowed. The right 
panels of Eigs. 2||4 show the extreme case where Cgjyj = 1. In 
this case the minimum possible radius is 9.0 km. 

Disconnected hybrid branches are of special interest, be¬ 
cause they give a characteristic signature in mass-radius mea¬ 
surements. Eor both the hadronic EoSs that we study, discon¬ 
nected hybrid star branches are excluded by the M^ax > 2 M 0 
constraint for Cgjy[ ^ 1/3, and even for larger Cgj^j they only 
arise if the hadronic matter EoS is stiff. Explorations of the ul¬ 
trastiff hadronic DD2-EV EoS indicate that disconnected hy¬ 
brid star branches can occur at moderate Cgj^, and it would be 
interesting to include this EoS in a future study. 

Our work is intended to motivate the use of the CSS 
parametrization as a framework in which the implications of 
observations of neutron stars for the high-density EoS can be 
expressed and discussed in a way that is reasonably indepen¬ 
dent of the modeling of the EoS of the high-density phase 
(quark matter in our case) |4| |26l . 

As an application to a specific model, we performed cal¬ 
culations for the ECM quark matter EoS. We showed that 
the ECM equation of state can be accurately represented by 
the CSS parametrization, and we displayed the mapping be¬ 
tween the ECM and CSS parameters. We found that ECM 
quark matter has a speed of sound in a narrow range around 
CgM = 0.3, and the ECM family of EoSs covers a limited re¬ 
gion of the space of all possible EoSs (Fig. [7|i. Once the ob¬ 
servational constraint Mmax > 2 M 0 is taken into account, the 
allowed region in the parameter space is drastically reduced 
to the shaded areas of Fig. This corresponds to the high- 
transition-pressure scenario, with a small connected branch of 
hybrid stars with tiny quark matter cores. Such stars would 
be hard to distinguish from hadronic stars via mass and radius 
measurements, but the quark matter core could be detectable 
via other signatures, such as cooling behavior. These hybrid 
stars have central densities larger than 6.5no in the BHF case, 
and 3.5no in the DBHF case (Fig.which means that, ac¬ 
cording to Refs. MMMM, hyperons could play an im¬ 
portant role in the BHF case, and they cannot be ruled out even 
in the relativistic DBHF case. As discussed in Sec. m we ig¬ 
nored hyperons because we are already using two different 
hadronic EoSs, one stiff and one soft, to estimate the sensi¬ 
tivity of our results to the hadronic EoS, and the effect of hy¬ 
perons remains unknown. In the future we hope that more ex¬ 
perimental data will constrain the high-density hadronic EoS, 
including the hyperonic contribution. 
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